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I A new analysis of S'-wave production amplitudes for the reaction 

' T^^Pi 'K~^'K~n on a transversely polarized target is performed. It is 

. based on the results obtained by the CERN-Cracow-Munich collab- 

\ oration in the vrvr energy range from 600 MeV to 1600 MeV at 17.2 

GeV/c 7r~ momentum. Energy-independent separation of the S-wave 
pseudoscalar amplitude (vr exchange) from the pseudovector amplitude 
r~| ! (oi exchange) is carried out using assumptions much weaker than those 

in all previous analyses. We show that, especially around 1000 MeV 
and around 1500 MeV, the ai exchange amplitude cannot be neglected. 
The scalar-isoscalar vrvr phase shifts are calculated using fairly weak 
\ assumptions. 

Below the KK threshold we find two solutions for the vr — vr phase 
shifts, for which the phases increase slower with the effective vr — vr 
mass than the P-wave phases. Both solutions are consistent with a 
broad /o(500) but only one is similar to the well-known "down" solution. 
We find also the third solution (with a somewhat puzzling behavior of 
inelasticity) which exhibits a narrow /o(750) claimed by Svec. All the 
solutions undergo a rapid change at the KK threshold. Above 1420 
MeV the phase shifts increase with energy faster than those obtained 
without the polarized-target data. This phase behavior as well as an 
increase of the modulus of the ai-exchange amplitude can be due to the 
presence of the /o(1500). 
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1 Introduction 



Study of scalar mesons is one of the central points of light quark spec- 
troscopy. In addition to ordinary qq mesons, some KK bound states [0, || and 
lowest -lying glueballs are expected as well ^ ^, 0. Unfortunately, the 



experimental situation is still far from being clear M,H, 10, 11|. In the effec 



tive mass region above 1000 MeV, a rich spectrum of scalar mesons has been 



recently proposed and discussed in many experimental |T^, |T3] and theoreti- 
cal 0, 1^, papers. Generally, the proposed scalar states have a coupling to 
the vrvr channel strong enough to manifest themselves through energy depen- 
dence of the TTTT interaction amplitudes. In the past few years scalar resonance 
/o(1500) related to a hypothetical lowest lying glueball state was announced 
H, ^ ITSf. Thus, a study of the vrvr interaction near 1500 MeV is important 
also in this context. Below 1000 MeV the status of scalar mesons is also un- 
clear since the existence of a broad cr or /o(750) meson still remains an open 



question |g, |16|, |17| |19 



One of the main sources of information on the production of scalar states 
is the TTTT partial wave analysis (PWA) yielding the S'-wave. It should be 
stressed that study of S'-wave objects does require the partial wave analysis 
to "subtract" contributions of leading mesons p(770) , /2(1270) and p3(1690) 
which dominate the total cross section. Virtually all PWA's were based on 
the old CERN-Munich experiment [pOl which supplied 3x10^ events of the 



reaction 

7T~P -K^-K^n (1) 

at 17.2 GeV/c. The number of observables provided by such experiment is 
much smaller than the number of real parameters needed to describe the partial 
waves. Consequently, the dominance of pseudoscalar exchange, equivalent to 
the absence of pseudovector exchange and several other physical assumptions 
have been made in previous studies p0|-p5|. 



In this paper we use results of PWA performed in the energy range from 
600 MeV to 1600 MeV (in 20 MeV bins) obtained with the help of the polarized 
target experiment. This experiment, performed 20 years ago by the CERN- 
Cracow-Munich collaboration, provided 1.2x10^ events of the reaction 

^ ~^ TT^TT^n (2) 

also at 17.2 GeV/c [Q. Combination of results of both experiments yields 



a number of observables sufficient for performing a quasi-complete and en- 
ergy independent PWA without any model assumptions. This analysis is only 
quasi-complete because of an unknown phase between two sets of transversity 
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amplitudes. Nevertheless, intensities of partial waves could be determined in 
a completely model-independent way. This removed ambiguities appearing in 
earlier studies, except for the old "up-down" ambiguity ||2l|. The "up" solu- 



tion contains an S'-wave resonance just under the p(770) and of similar width, 
while the " down" S'-wave modulus stays high and nearly constant all the way 
to the /o(980) . It was only Svec|TB, 17] who argued persistently in favour of 



the "up" solution, using both the 17.2 GeV data as well as data on the reac- 
tion vr+pi — > Tr~^iT~p at 5.85 and 11.98 GeV/c. However, general belief (see e.g. 
p7|] , [ p8[| ) was that the "up-down" ambiguity had been resolved definitely in 



favour of the "down" solution. We disagree with this belief since all the pre- 
vious studies of the full 17.2 GeV/c data were consistent with both the "up" 
and "down" solutions. The same is true for the present analysis. We stress 
this point because the mini-reviews in the last RPP editions (see e.g. Ref. 
p9|] ) contained the sentence "BECKER 79B jS^] excluded a narrow resonance 



behaviour for Sq ..." , contrary to what is stated explicitly in another paper of 
the same collaboration (see Sect. 6 of |2^): "...our results do not give a clear 
answer to this ambiguity...". 

In this paper we make another step in the analysis of 17.2 GeV/c data 
attempting to bridge two sets of transversity amplitudes. The phase of each 
S-wave transversity amplitude is fixed by requiring the phase of the leading P, 
D, F-waves to follow roughly the phase of the Breit-Wigner amplitudes of the 
p(770) , /2(1270) and p3(1690) resonances in the low, medium and high mass 
region respectively. This fairly reasonable assumption allows us to separate 
explicitly for the first time the pseudoscalar and pseudovector amplitudes in 
the S'-wave. 

In Sect. 2 we present mathematical formalism needed to separate the 
one-pion and ai exchange amplitudes for the reaction on a polarized target. 
Section 3 contains a short description of the PWA done by the CERN-Cracow- 
Munich collaboration for reaction (0). In Sect. 4 we present our analysis 
yielding pseudoscalar and pseudovector reaction amplitudes. Further on, from 
the pseudoscalar amplitude we extract the / = 0, S-wave amplitude describing 
the nil elastic scattering amplitude. Our results are discussed in Sect. 4 and 
summarized in Sect. 5. 

2 Amplitudes describing pion pair production 
in the / = 0, S'-wave 

a) Separation of pseudoscalar and pseudovector exchange amplitudes 

Let us denote by /o a system of two pions in a relative S-wave isospin 
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state. Transition amplitude for the /o production process vr p ^ fo n can 
be written as the following matrix element 

Ts^sr. =< <l^75 + ^B'y,^^{p^+pfr\u;i >, (3) 

where Pi,P2,Pn and pf are proton, neutron, incoming pion, and final /o four- 
momenta, Sp and Sn are proton and neutron spin projections, and 
are the corresponding four-spinors, A and B are functions of the Mandelstam 
variables s = (pi +Ptt)'^ and t = {pi —P2Y at fixed /o mass m-,tn- Part A of the 
amplitude corresponds to the pseudoscalar (or one pion) exchange while part 
B describes the pseudovector exchange which we shall briefiy call ai exchange 
since we expect that the ai meson exchange amplitude constitutes its major 
contribution. Functions A and B have to be determined from experiment. 
Using s-channel helicity amplitudes in the cm. 7r~p system one can derive the 
following two independent amplitudes ||2l[] : 



if++ = Tii =-T_i_i, (4) 
22 22 ^ ' 

H+- = Ti_i=T_ii, (5) 



22 22 



u I ^ r-, B m)-ml e, 

H++ = {7777V -tmin - 77 rZ7 ) COS — , (6) 

In (5-6) 0s is the neutron scattering angle (with respect to proton direction), 
M and m^r are nucleon (proton and neutron average) and pion masses, tmin 
and tmax are expressed by momenta pi, p^ and the corresponding cm. energies 
Ex, E2. 

trnin = 2{M^+PlP2-EiE2), (8) 
tmax = 2(M2 - p,p2 -E1E2). (9) 

The scattering angle is related to the four-momentum transfer squared t: 



Qg / t-min t , , 

^^^T = V t ■ -t ■ 

V ^mm ^max 

In this paper we use two amplitudes g and h, closely related to and 
adequate for describing /o production on a transversely polarized target: 

g=<ni \T\p t>= (ff+- - iH++) expi^Qg), (11) 

h=<n\ \T\p [>= {H+_ + iH++) exp{-Ues). (12) 
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In this case, nucleoli spin is quantized along the vector n = pi x p2 normal to 
the production plane (in (11) and (12), arrows t and | denote spin projections 
parallel or antiparallel to n). It can be shown that the remaining two matrix 
elements vanish: 

< n UT\p >=< n i \T\p i>= 0. (13) 

Separation of invariant amplitudes A and B for the reaction ir^p ^ fo n 
for fixed values of t and m.„T, can be done in the following way. Using (6,7) 
and (11,12) we define components qa, Qb and Ha, hs of amplitudes g and h as 
follows: 

A B 
9 = 9a + 9b, 9a = ^ U gB = —rV, (14) 

h = hA + hB, hA = ^U* hB = Y^*, (15) 



where 



U = {b-a)cs + i{bc'^ + as'^), (16) 
V ={l/a-l/b)cs-i{sya + cyb) (17) 



with the following shorthand notation: a = \/—tmaxi b = y/—tr, 



TTif — m 
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c = cosf^G,), s = sinf^Go 



From (14,15) one can obtain the desired amplitudes and ^r: 

^{gV*-hV), {li 



2M W 

^r = ^i-gU* + hU), (19) 

where 

W = UV* - U*V. (20) 

Solutions (18,19) are well suited for the analysis of experimental data taken 
with narrow t-bins. The CERN- Cracow-Munich data have been, however, 
averaged over a relatively wide t-range: 

h = -0.2 (GeV/c)2 <t<t2 = -0.005 {GeV/cf. (21) 

Therefore, we need amplitudes 'g and h averaged over this t-range, where 

g = -^ r git) dt (22) 



^2 ~ ^1 



5 



with similar equation for h. Then according to (14-15) we write 

g^J-AU+lrBV (23) 
1 



h = AU* + -rBV*, (24) 

2M 2 ^ ^ 



where symbols AU, BV etc. are defined as in (22). 

In order to proceed further we have to assume some t-dependence of func- 
tions A and B allowing for complete freedom in the effective mass dependence. 
Thus, we write 

A{mT,^, t) = ^o("T.7r7r)p("^7r7r, t) (25) 

and 

B{m^^, t) = Bo{m^^)q{m^^, t), (26) 

where p{mT^T^,t) and q{mT^-„,t) are postulated functions of t (which may also 
depend on nij^^), a4o(m7r7r) and i?o(T7^7r7r) depend on myrTr only and their be- 
haviour should be determined from experiment. A possible parametrization of 
p{m^^,t) is: 

p{m^^,t) ^ — 7, (27) 

— t 

which is equal to the pion propagator multiplied by the exponential form factor 
with parameter a{mT^T^) being an a priori unknown function of m^^. This shape 
of functional dependence has been introduced in many pion exchange models. 
For q{mTrTr,t) we have used four different parametrizations: 1, t, expipt) and 
texp{bt) with b—A GeV~^. As we shall show later, separation of averaged 
amplitudes g and h into a sum of averaged parts 'gX + 'gs and Ha + hs as 
in (14-15) and (23-24) is largely insensitive to the form of parametrization of 
function q{rnT,T,,t). Functions gAjhA, Wb and Hb are linear combinations of 
amplitudes g and h: 

gA^Cig + cji, (28) 
h2^dig + d2h, (29) 
gB = d2g- C2h, (30) 
hB^-dig + cih, (31) 

where the complex coefficients are 



ci = qV* pU/D, (32) 

C2 = -WpU/D, (33) 

di = qV^^/D, (34) 

d2 = -qVpU^/D (35) 
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and 

D = pUqV* -PU*W- (36) 

Having functions y and h experimentally determined, Eqs (28-36) make it pos- 
sible to calculate the unknown functions Aolrrij^T^) and BolrriT^^) in the following 
way: 



2M pU 

and 



(37) 



= w ^''^ 

From (25-26) we then obtain the desired functions A and B, so that separation 
of amplitude (3) into pseudoscalar and pseudovector parts is finally achieved. 

b) Determination of scalar-isoscalar pion-pion interaction ampli- 
tude 

Separation of the pseudoscalar amplitude, dominantly corresponding to 
the one pion exchange contribution to the n~p fo n process, is the first 
step in the determination of the scalar-isoscalar pion-pion amplitude Qq. This 
amplitude can be calculated if both the S'-wave tt+tt^ — > tt+tt^ amplitude as 
and the 1=2 S'-wave amplitude 02 are known: 

flo = 3 05 — -02 (39) 



(see |]21| p. 12 and notice that the /=!, S=0 contribution vanishes). Amplitude 
as is closely related to Aq/2M given by (37) since Aq is a factor responsible 
for the TTTT interaction in (25): 

= h==fWI7^ (40) 



where is the incoming vr" momentum in the n~p cm. frame, g,r is the 
final pion momentum in the /o decay frame, 5'^/47r = 14.6 is the pion-nucleon 
coupling constant, and / is the correction factor. In this factor the averaged 
t-dependence of the pion- nucleon vertex function and the off-shellness of the 
exchanged pion are included, which allows us to apply this formula to the 
analysis of the data taken in a wide t-region. 

When writing (^0]) we have assumed that absorption effects due to the 
final state interaction of the tttt with the outgoing neutron can be neglected. 
This assumption is supported by the results obtained in studies of absorption 
effects using the Regge phenomenology (|^). The point is that in the case of 
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dominance of the nucleon helicity flip amplitude, final state interaction effects 
are relatively weak. Thus, we can expect that errors caused by those effects 
are smaller than experimental errors. 

Amplitude 02 can be measured in the process tt^ p ^ 7r+7r+ n, provided 
the partial wave analysis is done and a similar separation of pseudoscalar and 
pseudovector contributions is performed. Although studies of tt^ti^ and tt^tt^ 
systems were in the past ||32|, the above-mentioned separation, which requires 
polarization measurements, has never been performed. Therefore, we have to 
rely on determination of the 1=2 amplitude based on the assumption that 
one-pion exchange dominates in the process under discussion. The 1=2 S- 
wave amplitude has been calculated using the data of |^ and the pion-pion 



separable potential model previously applied to the description of the coupled 
channel tttt and KK 1=0, S'-wave interactions [Q]. Here for the 1=2 channel 
we use a very simple two-parameter pion-pion potential of rank one: 

V^{p,p') = XG{p)G{p'), (41) 

where A is a constant and 



Mtt 1 

^(^) = V ^T^ (42) 

V p + p 

is a form factor with one range parameter f3\ p and p' are pion cm. mo- 
menta in the initial and final states respectively. The calculated phase shifts 
82 corresponding to the amplitude 

02 = ?>m.52exp{i52) (43) 

for A = = -0.1309 and /3 = 3.384 GeV are shown in Fig. 1. The analytical 
form of the amplitude 02 can be found in Appendix A of P]. This amplitude, 
along with as obtained from (40), was used to calculate oq given by (39). 
Amplitude 02, being generally smaller than 05, cannot be, however, neglected 
in (39) as shown in Sect. 4 d. Amplitude is normalized to Argand's form: 

,44) 

where 5 is the /=0, S=0 phase shift and r] is the inelasticity coefficient. 



Model-independent determination of partial- 
wave amplitudes 

This section is a short recapitulation of what was extensively described in 



the old papers of the CERN- Cracow-Munich collaboration [^, to 
which the reader is referred for more details. 
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In the analysis, 3 x 10^ events on a hydrogen target ||20(] were combined with 
1.2 X 10^ events on a polarized (butanol) target, both for 17.2 GeV/c vr". The 
former events yield moments while the latter provide the p\j and r^^ mo- 
ments from interactions with protons of hydrogen bound in butanol molecule 
(protons in carbon and oxygen nuclei cannot be polarized). In this way we 
obtain a quasi-complete description of the reaction: 

7r"p| n^Ti^n. (45) 

In the case of a transversely polarized target the differential cross section can 
be written as 



dmT,T,dtdVt ^^j^ j^^^ j^j^ 

(46) 

where: 

m^^ - effective mass of the vr+vr" system, 

t - four-momentum transfer squared from the vr" beam to the tt^tt" system, 

- spherical harmonics, 
fl - decay angles of the tt" in the t-channel 7r"'"7r~ rest frame, 
ijj - polarization angle, 
P - degree of polarization. 

Table 1 



Number of parameters and observables in partial-wave analysis, / denotes the vr 
orbital momentum in the vr+vr" rest frame and m ~ its projection. 



^max(wave) 


1(P) 


2iD) 


3(F) 


Number of amplitudes (m < 1) 


8 


14 


20 


Number of real parameters (m < 1) 


14 


26 


38 


Number of t^f moments (M < 2) 


6 


12 


18 


Number of moments (M < 2) 


6 


12 


18 


Number of moments (M < 2) 


3 


7 


11 


Total number of moments (M < 2) 


15 


31 


47 



Moments of angular distribution t^, and are bilinear combina- 
tions of partial wave amplitudes. As described in the previous Section, we use 
nucleon transversity amplitudes g and h corresponding to a given naturality 
exchange. It should be stresssed that the number of moments is smaller (see 
Table 1) than the number of real parameters characterizing the amplitudes. 
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and therefore model-independent partial-wave analysis is not possible. Con- 
sequently, all TT+vr" phase shift studies that do not use polarized target data 
were based on non-trivial physical assumptions like vanishing of spin-nonfiip 
amplitudes (in our language this corresponds to gi = hi, where i denotes l,m) 
in the unnnatural spin-parity amplitudes and phase coherence between the 
m = and m = 1 amplitudes. It was shown in |TH| , ^ that these assump- 
tions are badly broken by the polarized-target data. Unfortunately, this fact 
has been ignored in all subsequent vrvr studies with the notable exception of 



Svec papers [17, IS 



On the other hand, additional knowledge of pj^ and rj^ moments yields 
the total number of observables which, as seen in Table 1, slightly exceeds the 
number of real parameters. Since the present analysis is restricted to low t, we 
can ignore all m > 1 amplitudes as all M > 2 moments vanish (for a high-t 
study see [p5| ). 

Terms appearing in amplitude combinations are of the type \gi\'^, 
Re{gig*) or Re{hih*) but there is no mixed term like Re{gih*). Consequently, 
we cannot determine the relative phase between the g and h amplitudes in 
a model-independent way but the transversity amplitudes can be completely 
determined independently within each set. This includes moduli and relative 
phases with the warning that all the relative phases within each set can be 
multiplied by —1. 

In the analysis we expect to find some discreet ambiguities that arise from 
bilinearity of the equations. In order not to lose any solution, great effort was 
devoted to providing many different starting points to the MINUIT program 
[p^] . They were as follows (see |Q for details): 

i) exact analytical solutions (possible for Imax = 1 only), 

ii) approximate analytical solution assuming that intensity of the Pq wave is 
much smaller than that of the 5* and Do^waves (for l^ax = 2), 

iii) approximate analytical solution assuming phase coherence (for l^ax = 2 
and Imax = 3), 

iv) several sets of random values, 

v) results of the fit from neighbouring bins. 

It should be stressed that approximations in ii) and iii) were used for finding 
the starting values only; it was always the exact formulae that were fitted. In 
[EB] such fits were done in 40 MeV bins for 



580 MeV < m^^ < 1780 MeV 
0.01 GeVVc^ < \t\ < 0.20 GeVVc^. 

Later, similar analysis was performed in finer mass bins (A m^rTr = 20 MeV) 
in a slightly different kinematical region, i.e. 

10 



600 MeV < m^^ < 1600 MeV 
0.005 GeVVc^ < \t\ < 0.200 GeVVc^. 

The fits became more difficult and time-consuming as higher partial waves 
were included; therefore, the MINOS error analysis (see for details) was 



required for all parameters for l^ax = 1 only. For higher masses, when more 
waves were needed, the MINOS errors were calculated for the leading m = 
waves only. The solution was considered to be acceptable only if the MINOS 
error analysis was possible; the values, being generally good, were hardly 
helpful in selecting solutions. Quite often there was only one such solution in 
a given bin, although many different starting points were used. However, in 
the mass region below /o(980) there are two branches of solutions, best seen 
in the moduli of the S-wave g and h amplitudes around 900 MeV. This reflects 
the old "up-down" ambiguity. 

Intensities li = \gi\'^ + \hi\'^ of partial waves obtained in the 20 MeV analysis 
were published in in fact they still represent the most accurate measure- 
ments of the /2(1270) parameters Moduli and relative phases of transver- 



sity amplitudes are used in this paper for the ffist time. 



4 Results 

a) Determination of the S— wave 

In Figs 2-4 we show the experimental results obtained by the CERN- 
Cracow-Munich collaboration as described in the previous Section. Indepen- 
dent variables used in their partial wave analysis were: the sum + 
(Fig. 2a), the ratio (Fig. 2b), as well as the phase differences 'd^ — -d^ 

and — 'i?^ (L, L' = S,P,D,F) (Figs 3 and 4). In our analysis we have 
assumed that phases of the partial waves are described mainly by phases of 
the final state interaction amplitude of the vrvr system. It means in particular 
that phases of the P, D and F-waves follow phases of p(770) , /2(1270) and 
P3(1690) decay amplitudes into tttt . Our PWA has been done in three m^^r 
effective mass regions from 600 MeV to 1600 MeV. 

i) 600 MeV - 980 MeV 

In this region one can safely assume that only the S and P-waves contribute 
since the D-wave is weak even at the upper limit of this region. It is only in 
this region that fully analytical solutions of the PWA equations are possible 
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0^] . The PWA analysis however, yields two solutions ("up" and "down"' 



which are distinctly different in the m^^r region from 800 MeV to 980 MeV. In 
the "up" solution the sum | ^ p + | /i p exhibits a maximum for mTTTr ~ 770 
MeV, but in the "down" solution the moduli of | ^ p and \ h are roughly 
constant from about 750 MeV to 980 MeV. The shape of the S'-wave "up" 
modulus was used by Svec to claim the existence of the /o(750) meson [|T^ |T7|. 



S'-wave phases of reaction have been determined from phase differences 
between S and P-waves. In addition to the "up-down" ambiguity in the 
moduli of the g and h transversity amphtudes, there is also a phase ambiguity 
in each mj^T^ bin. This ambiguity comes from the mathematical structure of the 
PWA equations from which only cosines of the relative phases of the partial 
waves of reaction (Q) can be obtained. In our analysis we present two arbitrary 
choices of the S'-wave phases which form data sets shown in Figs 3 and 4. In 
the first set, called "steep", S'-wave phases grow faster than P-wave phases. 
In the other set, called "fiat", increase in S'-wave phases is slower than that for 
P-waves. Thus two sets of possible phases ("fiat" or "down") combined with 
two branches of moduli ("up" and "down") lead us to consideration of four 
solutions for the amplitudes g and h which can be called "up-steep" ," down- 
steep", "up-fiat" and "down- fiat". Following this splitting of the amplitudes 
also the vr — tt phase shifts which will be derived from them will be similarly 
labelled. In order to avoid a possible future confusion let us remark here 
that in many papers in past the words "up" and "down" served mainly to 
distinguish the S-wave phase-shifts: the values of the phase shifts "up" were 
larger than the corresponding values of the solution " down" . However, in the 
models in which the ai exchange has been neglected and the S-wave amplitude 
was assumed elastic, the modulus of the "up" amplitude vr^p — > vr^Tr+n was 
sma//er than the modulus of the amplitude "down" for the effective t[~tt^ mass 
larger than about 800 MeV. The reason of this correlation can be understood 
if we remember that the phase shifts of the " down" solution were close to 90° 
while the phase shifts of the "up" solution were considerably larger in that 
effective mass range and the modulus of the amplitude was proportional to 
I sin 5 |. In our case, however, we do not neglect an important contribution of 
the ai exchange, so we have more combinations of the possible solutions for 
the amplitudes g and h. 

ii) 980 MeV - I46O MeV 

Here, the S, P and D-waves have been included in the PWA analysis. In 
this region the /2(1270) production is very strong so the D-wave dominates 
and the S'-wave phases have been evaluated from phase differences between 
the S and D-waves. In our model we have assumed that in the 980 MeV - 
1460 MeV region differences between the P and D-wave are positive since the 



12 



phase of the p(770) meson decay amphtude into two pions is larger than the 
corresponding phase of the /2(1270) meson amphtude. This allowed us to fix 
the sign of phase difference between the S and D-waves in each energy bin, 
thus avoiding phase ambiguity present in the first region. 

tu) 1460 MeV - 1600 MeV 

Here, the S, P, D and F (dominated by the p3(1690) meson) waves have 
been included in the PWA analysis. Consequently, the 5'-wave phases have 
been calculated from phase differences between the S and F-waves. The sign of 
each difference has been fixed by the assumption that phase differences between 
the P and F-waves as well as between the D and F-waves are positive. 

In each of the three regions production amplitudes of the p(770) , /2(1270) 
and Pa (1690) resonances follow the Breit-Wigner parametrisation of [^, |3^ , 
namely 

M = VA^ ^^'""" Sl-, i = 1, 2, 3, (47) 



where 




mJ Di{qr) 
and the Blatt-Weisskopf functions Di{qr) have the form: 



(4J 



Di{qr) = 1 + {qrY for P-wave, 

D2{qr) = 9 + 3(gr)^ + (gr)^ for D-wave, 

D^iqr) = 225 + 45(gr)2 + 6(gr)^ + (gr)^ for F- wave. (49) 

In (|47|) A is a normalization constant , m^, Tji and xpi are mass, width and 
inelasticity of the resonance R (-R=p(770) , /2(1270) or p3(1690) ), q is mo- 
mentum of any pion in the vrvr rest system, and qn stands for momentum q for 
m^TT = f^R- Range parameter r for p(770) is equal to 4.8 GeV~^ for the "up" 
solution, and 5.3 GeV~^ for the "down" solution |3^. For the /2(1270) and 



P3(1690) resonances r equals 10 GeV~^ and 3 GeV~^ respectively 

Simultaneous presence of the P and D-waves for m^yr > 980 MeV allowed 
us to check the validity of parametrisation of the P and D-wave phases by 
resonant amplitudes ([471) . In Fig. 5 we have shown phase differences between 
the P and D waves for 'g and h transversity amplitudes. The P and D-waves 
in the h transversity amplitude are well described by the p(770) and /2(1270) 
resonant amplitudes but in the case of the g amplitude such parametrisation 
is not sufficient. 

Since in the PWA analysis there is no phase relation between the g and h 
transversity amplitudes, we have defined the S'-wave phases in the following 
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way: 



= <! < - + 9M1270) 
^9^-< + ^P3(1690) 



for 600 MeV < m^^ < 980 MeV, 
A for 980 MeV < m^^ < 1460 MeV, 
for 1460 MeV < m^^ < 1600 MeV, 



(50) 



+ ^p(77o) + A for 600 MeV < m^^ < 980 MeV, 
^1-^^ + ^2(1270) + A for 980 MeV < m^^ < 1460 MeV, (51) 
- < + ^p3(i69o) for 1460 MeV < m^^ < 1600 MeV 



In (|50|) and ([5T|) Op^uQ), 6'/2(i27o) and 6'p3(i69o) are phases of the resonant amph- 
tudes defined in (^). Function A has the form: 

r 50.37° for 600 MeV < m^^ < 990 MeV, 

A = I 50.37° - 0.116°(m^^ - 990 MeV) for 990 MeV < m^^ < 1420 MeV, 
[ for m^^ > 1420 MeV. 

(52) 

Function A has been introduced in order to parametrize the differences between 
(i)^ — d^) and {dg — ^^)- We have determined this function empirically from 
a linear fit in the range of 980 MeV < m^^ < 1300 MeV. 

b) Separation of the 5— wave into pseudoscalar- and pseudovector- 
exchange components 

Transversity amplitudes 'g and h have been averaged over t as it was 
described in Sect. ||. Coefficients ci,C2,(ii and d2 depend on averages pU and 
qV , and therefore they also depend on the form of functions p{t) and g(t). We 
have checked the dependence of amplitudes (a = A, B) and A, B on 

the value of parameter a (Eq. for < a < 4 GeV~^. For both "up" and 
"down" solutions and for both "fiat" and "steep" sets of the S-wave phases 
below 1000 MeV, the differences in moduli and phases of amplitudes Ha 
and 'qbi are not higher than 3% for 600 MeV < m^^r < 1600 MeV. In the 
case of the A amplitude the differences between different solutions can reach 
10%. We have evaluated parameter a to be 3.5 GeV~^ from experimental data 
pO|] in the mass range of 710 MeV < m^^ < 830 MeV. Since we do not know 
its mass dependence outside this region, we have kept it fixed to 3.5 GeV~^ in 
the whole mj^^ range. 

We have also checked the dependence of pseudoscalar and pseudovector 
exchange amplitudes on the form of function q{t) for the following functions: 
q(t) = l,t, e^* and q(t) = te''*, where 6 = 4 GeV~^. The second and fourth 
forms come from parametrization of the Regge propagator for the ai exchange 
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at t ~ fST]]. Differences in moduli and phases of the A, 'gj, Ha-, 'Qb and Kb 
amphtudes caused by different functional forms of q{t) are smaller than 5%. 
Very small changes of amplitude B (smaller than 0.5%) are due to the strong 
factorization of the average qV ~ qV . Therefore, in the numerical calculations 
we have chosen q{t) = 1. 

Parameter / introduced in (^) has an influence on amplitudes and as 
defined in Sect. 0. We have evaluated this parameter by minimizing differences 
1 - r/ for 600 MeV < m^^ < 990 MeV for each solution ("up" and "down") 
combined with the phase set ("fiat" and "steep") separately. 

Since the dependence of transversity amplitudes and ha (a = A, B) 
on the shape of functions p{t) and q{t) is weak, we have used them to separate 
contributions of one-pion and ai exchanges. Amplitudes and depend on 
moduli and relative phases of amplitudes ^ and h . Since for m^rTr < 980 MeV 
there are two solutions for moduli ("up" and "down") and two solutions for 
relative phases ("flat" and "steep"), we obtain four combinations of amplitudes 
which will be further labelled "up-fiat" , "up-steep" , "down-fiat" and "down- 
steep" . In Fig. 6 and 7 the moduli and phases of these four amplitudes are 
shown. 

c) Properties of pseudoscalar- and pseudovector-exchange ampli- 
tudes 

Separation of amplitudes leads to the equality of corresponding moduli: 
I hA 1 = 1 gA I and | hs | = | 'gs \ (see Eqs |1^ and |T5[). However, phases of Ha 
and Kb differ from the corresponding phases of and . Their differences, 
being functions of the effective mass, are related to the phases of complex 
coefficients Cj and rfj (i = 1, 2 in Eqs p8|-^). The first phase difference, defined 
as the phase of Ha minus the phase of 'gX, increases monotonically from about 
10° to 50° in the effective mass range from 600 MeV to 1600 MeV. The second 
phase difference, defined as the phase of hs minus the phase of 'gs, is almost 
constant (~ —174°) in the whole region of the effective mass. This behaviour 
follows from the fact that the imaginary part of coefficient U (Eq. 0) is 
much smaller than the real part of U, and the imaginary part of coefficient 
V dominates over its real part (Eq. |1^. Since we assume that functions p(t) 
and q{t) are real, the phase of amplitude A is the arithmetic average of the 
gy^ and Ha phases while that of amplitude B is the arithmetic average of 
the g^ and hs phases (compare Eqs|l^ and ^). Their values differ strongly 
between four sets of solutions. The g^ phases systematically increase for 
the nil effective mass between 600 MeV and about 980 MeV. Sudden phase 
change and strong reduction of the 'g^ moduli for m^rTr — 1000 MeV results 
from the appearance of the narrow resonance /o(980) . Similarly decrease in 
the phases and moduli of amplitudes for ninn > 1400 MeV can be related 
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to the presence of another scalar resonance of mass between 1400 MeV and 
1500 MeV. In most cases, errors of the phases and moduh of are smaller 
than the corresponding errors for . 

The Oi exchange amplitude 'g^ amounts on average to about 20% of the 
pion exchange amplitude but around 1000 MeV and 1500 MeV it is of the 
same order as . A slight increase in the modulus and a decrease 
in the g^^ modulus for the effective mass above 1420 MeV (seen in Figs 6 
a,b and 7 a,b) can be related to the presence of a scalar resonance /o(1500) 
coupled to channels such as nn, rjrj, rjrj' and in (or aa) 0, |]10| - [1^. On 



the basis of the pp 57r° annihilation data, S. Resag from the Crystal Barrel 
Collaboration [|lT] suggested a possible coupling of the /o(1500) to 7r(1300)7r, 
where 7r(1300) couples in turn to (vrvr)^?? with (TrTr)^ denoting the vrvr pair in 
the / = 0, S'-state. The possible enhancement of the | 'gs \ amplitude seen in 
Figs 6 and 7 above 1400 MeV can be, however, attributed to another decay 
mode of /o(1500) , namely to a system aivr, where ai — » (7r7r)57r. However, the 
small value of the partial decay width of ai into (TrTr)^'/?, which - according 
to the Particle Data Group |^ - is equal to 0.3% of the total ai width, is 
questionable. This number follows from just one analysis related to the 



not well confirmed four-quark model of Jaffe The data gathered by the 



ACCMOR collaboration ||39[ for reaction n^p 37rp indicate a much larger 
ai coupling to (7r7r)57r (even several tens in percent of the total ai width). 
In reaction n'p — > n^TT^n a possible subprocess is the Oi exchange followed 
by interaction 7r~ai {7t^tt~)s in the isoscalar state. It seems therefore the 
pseudovector exchange process can also contribute to possible production of 
the /o(1500) scalar resonance in the reaction under consideration. 

d) Pseudoscalar / = 0, S— wave amplitudes 

Finally, let us extract the vrvr scalar-isoscalar amplitude Oq from pseu- 
doscalar amplitude g^ just separated from pseudovector amplitude g^ . We use 
formulae ( p7| ) and (p9D-(^) to obtain phase shifts 6 and inelasticity coefficient 
7] corresponding to Oq. In Fig. 8 we show the dependence of those parameters 
on the effective mass for solutions labelled "steep" (corresponding to a fairly 
steep behaviour of the phases of amplitudes g and h ). For both solutions 
("up-steep" and "down-steep" ) we observe a characteristic fast increase of 
phase shifts near 760 MeV (see Figs 8 a,b). Now all the solutions have to pass 
a test connected with the fitting of the complex coefficient / =| / | e*"^ to 
the requirement that the inelasticity coefficient rj = 1 in the whole mass range 
from rriTrn = 600 MeV to the KK production threshold. 

It appears immediately that the the "down-steep" solution fails this test. 
Such an attempt leads to unnacceptably high /NDF = 80/17 value ob- 
tained in the best case for | / |= 0.78 and ip = —8.4°. In addition, there 



16 



is a strange behavior of rj i.e. a clear violation of unitarity constraint below 
m,r7r < 720 MeV and above 820 MeV while rj < 1 between these two values 
(see Fig. 8 c). It should be noted that if we changed normalization to keep 
f] < 1 in the whole mass range in question we would deal with a puzzling large 
inelasticity between m^^ = 720 MeV and 820 MeV. Therefore, in our opinion 
the "down-steep" solution is nonphysical and should be excluded. 

The situation is different for the "up-steep" solution, for which we can find 
a constant factor such that rj does not deviate from unity {x^/NDF = 15/17) 
for I / 1= 0.73 and (p = —9.2°. However, also here we observe a puzzling be- 
haviour with f] systematically exceeding unity for mT^Tr < 720 MeV and above 
m^^ = 820 MeV but staying below 1 in the intermediate mass range. 

On the other hand both "flat" solutions exhibit natural behavior with rj 
fluctuating around unity. In addition, minimizing the t] — 1 differences be- 
low the KK production threshold we obtain reasonable values of / close to 
1, namely | / |= 0.89, if = -4.4° for the solution "up" and | / |= 0.84, 
ip = —17.8° for the solution "down". Thus we are left with two favoured, nicely 
behaving, solutions ("up-fiat" and "down- fiat") and one, somewhat queer, but 
still acceptable "up-steep" solution. Let us remark that on/?/ the last solution is 
consistent with the relatively narrow /o(750) meson Q claimed by Svec |0 
on the basis of moduli of transversity amplitudes determined with the help of 
the polarized-target data. 

In all solutions there is a sudden drop in rj for the effective mass near 1000 
MeV, caused by an opening of a new KK channel. Another decrease of t] can 
be seen above 1500 MeV. 

In Figs. 8 a,b and 9 a,b the phase shifts S for all solutions have been com- 
pared with those obtained from the analysis of the n^p Tr+vr^n reaction 
on an unpolarized target [20] (solution B), where separation of the vr-exchange 
and ai-exchange amplitudes was impossible. In Fig. 9 b in the mass region 
from 600 MeV to about 1400 MeV we see only minor differences between phase 
shifts corresponding to the "down- fiat" solution and the results of [20]. In the 
"up-flat" solution around 900 MeV, however, the values of 6 are higher than 
those of ref. [20] by several tens of degrees as seen in Fig. 9 a (a qualitatively 
similar difference led many years ago to the "up-down" nomenclature). The 
"up-steep" solution (see Fig. 8 a) with a narrow resonance is entirely different 
from the results of ref. [20]. Let us notice that the characteristic "jump" of the 
"up-steep" phase shifts at 980 MeV is smaller (about 100°) than those in the 
the "up-fiat" (^ 120°) and "down-fiat" (^ 140°) solutions. 

Finally let us comment on the behaviour of the phase shifts at higher 
effective masses, well above the i^'i^'threshold. Here the different solutions 

^Following the PDG convention we call this object /o(750) instead of (7(750) used in the 
original papers. 
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are quite close each other. This follows from our assumption about the phase 
difference between the D- and S- waves (formulated in Sect. 4a ii) and rather 
small differences between the phases ^ of the factor / (eq. 40). 

A short comment on the comparison of errors obtained for the phase shifts 
in this study and errors obtained in is in order here. Errors in the partial 



wave analysis have been obtained from the MINOS subroutine of the 
MINUIT program used to fit data independently in each effective mass bin. 
Although the errors of the B analysis in |^ are smaller, they do not include 
systematic effects corresponding to simplifying assumptions like, for example, 
dominance of the s-channel nucleon helicity flip amplitude and phase coherence 
of unnatural spin-parity exchange amplitudes with helicity and 1. In our 
analysis we avoid these assumptions which have been found to be badly broken 



by the polarized target data already in [p^ j , pO| , [p^ . 

Finally, let us note that the / = 2, S'-wave amplitude 02 plays an important 
role in the determination of amplitude ao (Eq. In particular, especially at 
high effective masses, it influences the inelasticity parameter 77 very much. As 
described above, phases of amplitude 02 have been obtained in the experiment 
on an unpolarized target . Thus, they contain some unknown errors related 
to the unknown contribution of the a\ exchange. We feel, however, that those 
errors can be less important than other errors of amplitude as calculated from 
gj, (EqsHandP). 

e) Resonance interpretation of results 

Let us now discuss behaviour of phase shifts and inelasticity coefficients 
in terms of the scalar /=0 resonances, both already known and those newly 
postulated (see Figs 8 and 9). 

Systematic increase in phase shifts below 1000 MeV can be related to the 
existence of a broad o meson which we have found in and called /o(500). 
In we analysed data obtained on an unpolarized target. Applying now 
the same model to the "down-fiat" and "up-fiat" solutions obtained for the 
polarized target we obtain parameters of the /o(500) meson very similar to 
those given in [Q. In this way, presence of /o(500) is reinforced. We should 
mention here that N. A. Tornqvist and M. Roos [0 also support the existence 
of the a meson naming it tentatively /o(400 — 900). The a meson of mass 555 
MeV has also been reported recently by Ishida et al. ([|31) analysis of 

the TTTT phase shifts especially near the pion-pion threshold. The threshold 
behaviour of the scalar /=0 scattering amplitudes is very much influenced by 
the presence of the /o(500) state (see |^). The a meson plays also a very 
important role in nuclear physics and in the description of nucleon-nucleon 
interactions. For example, in |^ typical value of 550 MeV has been used for 



its mass, which is in agreement with the above findings. 
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Since our "up-steep" solution (see Sec. 5d) corresponds to the narrow 
/o(750) claimed by Svec on the basis of (mainly) the same data let us comment 
on the similarities and differences of two approaches. Svec uses an " analytical 
solution+Monte Carlo" method to extract amplitudes from experimental mo- 
ments while our amplitudes come from the standard minimization method. 
The advantage of the latter is an applicability to higher mass range where 
more partial waves contribute and analytical solution is not possible. In the 
very recent draft he compares both methods for 171,^ < 900 MeV finding 
an excellent agreement between two methods. But there is an important dif- 
ference in the further analysis. Svec finds narrow /o(750) in only one of two 
transversity amplitudes and argues against summing the squares of two ampli- 
tude moduli, which - according to him - distorts or even hides the /o(750). In 
our analysis we extract the dominating vr-exchange amplitude and weaker but 
non-negligible ai-exchange amplitude from connecting both transversity ampli- 
tudes. Studying the pure 7r-exchange amplitude we find one, less favoured but 
still possible, solution exhibiting the narrow /o(750). A simple Breit-Wigner 
fit with the energy-independent background yields m = (770 ± 8) MeV and 
r = (120±9) MeV, well consistent with m = (753±19) MeV and T = (108±53) 
MeV in Ref. UTSf. This solution, however, cannot be well described by the 
above-mentioned model of Ref. [2] without a substantial modification of the 
pion-pion interaction. 

As can be seen in Fig. 9 c, a sudden drop in inelasticity for the effective 
mass near 1000 MeV is caused by the opening of a new KK channel. This 
fact, along with a characteristic jump of phase shift 6 coupled with a rapid 
decrease in elasticity 77 near 1000 MeV (Fig. 9 a,b), is due to the narrow 
/o(980) resonance [0, We stress that this jump, although smaller than 

in two "flat" solutions, is seen also in the "up-steep" solution, since in some 
older studies the "up" solution was rejected as inconsistent with the narrow 
/o(980) . Analysing the polarized target data Svec has shown that an intensity 
of one S-wave transversity amplitude can be simultaneously parametrized in 
terms of two interfering resonances /o(750) and /o(980) |18| |. 

Another decrease in rj near 1500 MeV - 1600 MeV can be related to the 
opening of further channels like Air ( a a or pp ), rir]' or uu. For 7717^77 larger 
than 1470 MeV, phase shifts for both solutions ("up-fiat" and "down-fiat" ) 
show a systematically steeper increase than phase shifts corresponding to data 
obtained on the unpolarized target. Both facts may be related to the possible 
existence of the /o(1500) resonance 0, |^, [0-[jl2 . 
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5 Summary 



A formalism permitting extraction and separation of the .S-wave pseu- 
doscalar and pseudovector exchange amphtudes in the reaction vr^pj tt'^tt'ti 
on a transversely polarized target has been presented. A new analysis of the 
CERN-Cracow-Munich collaboration data obtained with 20 MeV bins at 17.2 
GeV/c 7r~ beam momentum on the polarized target has been performed in 
the TTTT energy range from 600 MeV to 1600 MeV. Already at the level of 
moduh this analysis yields two solutions between 600 MeV and 980 MeV (the 
"up-down" ambiguity) . 

Due to the lack of information on the relative phases between transversity 
amplitudes 'g and h we have assumed that the relative phases of both S*- 
wave transversity amplitudes g and h are generally governed by the phase 
behaviour of the dominant resonant P, D and F partial wave amplitudes 
corresponding to the p(770), /2(1270) and p3(1690) resonances decaying to the 
TTTT pairs and interfering with the S'-wavc. This leads to an additional twofold 
ambiguity since the relative phases can be cither added or subtracted. Thus we 
have "down- flat", "down-steep", "up-flat" and "up-steep" solutions. However 
the "down-steep" solution is shown to violate unitarity and can be ignored. 
The remaining three solutions are acceptable although the the " up-steep" one 
exhibits a peculiar behavior of inelasticity. 

The Oi exchange amplitudes are especially important at 1000 MeV and 1500 
MeV and cannot be neglected with respect to the tt exchange amplitudes. This 
puts in serious doubt all the PWA results which assumed absence of the ai 
exchange. 

Separation of the 7r-exchange from the ai-exchange allowed us to calculate 
the / = 0, S'-wave tttt amplitudes in a weakly-model-dcpcndcnt manner. In 
the low-mass region both "flat" solutions arc consistent with /o(500) while 
the phase behavior of the "up-steep" one agrees with the narrow /o(750). Up 
to the energy of about 1420 MeV, phase shifts of our "down-flat" solution 
agree within the errors with the one obtained without the polarized-target 
data. However, above 1420 MeV the phase shifts in all solutions increase with 
energy faster than those obtained without the polarized-target data. This 
phase behavior as well as an increase of the ai-exchange amplitude can be due 
to the presence of the /o(1500). 
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FIGURE CAPTIONS 



Fig. 1. The 1 = 2, S-wave tttt phase shifts versus the effective tttt mass. 
The curve represents the fit to the data [32]. 

Fig. 2. Results of the partial wave analysis of the CERN-Cracow-Munich 
collaboration: a) | ^ P + \h\'^,h) ratio | ^ | / | 7i |. Below 980 MeV, full 
and open circles represent the "down" and "up" solution, respectively. 

Fig. 3. Phase differences of the g transversity amplitudes obtained in 
the partial wave analysis of the CERN-Cracow-Munich collaboration: 
a) Solution "up" : phase differences ■j?^ — ■j?^ between the S and P-waves 
for the "flat" set (full circles) and the "steep" set (open circles), phase 
differences 'dg — between the S and D-waves (diamonds) and phase 
differences i?^ — i9g between the S and F-waves (squares), b) Solution 
" down" : notation as in a) . 

Fig. 4. Phase differences of the h transversity amplitudes obtained in 
the partial wave analysis of the CERN-Cracow-Munich collaboration. 
Notation as in Fig. 3. a) Solution "up", b) Solution "down". 

Fig. 5. Phase differences between the P and D-waves for the g {d^ — 
•d^ : open circles) and h {d^ — i?^ : full circles) transversity amplitudes 
versus the effective tttt mass. Dashed line represents the effective tttt 
mass dependence of function A obtained from a fit to differences (i?^ — 
d^) — {dg — d^) denoted by triangles. Solid line represents differences 
between phases of the p(770) and /2(1270) decay amplitudes. 

Fig. 6. a) Moduli of pseudoscalar | ^ | = | /i^ | (open circles) and pseu- 
dovector | ^ | = | /?._b | (diamonds) exchange amplitudes as a function of 
the effective tttt mass for the "up-fiat" solution, b) Same as in a) for 
the "up-steep" solution, c) Phases of pseudoscalar exchange amplitudes 
gj^ (open circles) and pseudovector amplitudes 'gs (diamonds) versus the 
effective tttt mass for the "up-flat" solution, d) Same as in c) for the 
"up-steep" solution, e) The phases of the pseudoscalar exchange ampli- 
tudes (open circles) and the pseudovector amplitudes hs (diamonds) 
versus the effective tttt mass for the "up-flat" solution, f) Same as in e) 
for the "up-steep" solution. 

Fig. 7. a) Moduli of pseudoscalar | ^ | = | Ha \ (full circles) and pseu- 
dovector I ^ 1 = 1 /is I (diamonds) exchange amplitudes as a function 
of the effective tttt mass for the "down-flat" solution, b) Same as in 
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a) for the "down-steep" solution, c) Phases of pseudoscalar exchange 
ampUtudes 'gj (full circles) and pseudovector amplitudes qb (diamonds) 
versus the effective tttt mass for the "down-fiat" solution, d) Same as 
in c) for the " down-steep" solution, e) Phases of pseudoscalar exchange 
amplitudes Ha (full circles) and pseudovector amplitudes Kb (diamonds) 
versus the effective tttt mass for the "down-flat" solution, f) Same as in 
e) for the "down-steep" solution. 

• Fig. 8. a) Scalar-isoscalar tttt phase shifts Sq as a function of the effective 
nn mass for the "up-steep" solution (open circles) and for data [20] 
(triangles), b) Same as in a) for the "down-steep" solution (full circles), 
c) Scalar-isoscalar tttt inelasticity coefficient r] versus the effective tttt 
mass for the "down-steep" (full circles) and "up-steep" (open circles) 
solutions. 

• Fig. 9. a) Scalar-isoscalar tttt phase shifts Sq as a function of the effective 
TTTT mass for the "up-flat" solution (open circles) and for data [20] (tri- 
angles), b) Same as in a) for the "down-flat" solution (full circles), c) 
Scalar-isoscalar tttt inelasticity coefficient rj versus the effective tttt mass 
for the "down-flat" (full circles) and "up-flat" (open circles) solutions. 



25 



-40 



400 600 800 1000 1200 1400 1600 

mTTTT (MeV) 



Fig.l 




\h\ 



2.5 
2 

1.5 
1 

0.5 




b) 




1 



7.1 



5.6 



4.9 



3.6 



800 1000 1200 1400 1600 

mTTTT (MeV) 



Fig. 2 




Fig. 3 



bJO 



200 
150 
100 
50 



a) 



<> 



i4 



-50 
-100 



4m 







<> 



CP 







><> 



<> 



800 



mTTTT (MeV) 



bJO 



200 
150 
100 
50 



b) 



Q 



-50 
-100 







o 



Q T 



<>< 







<> 




V 



><> 











CP 



1000 1200 1400 1600 



[p 







800 1000 1200 1400 1600 

mTTTT (MeV) 



Fig. 4 




Fig. 5 




800 1000 1200 1400 1600 

mTTTT (MeV) 




800 1000 1200 1400 1600 

mTTTT (MeV) 



Fig. 6ab 



bJO 




800 1000 1200 1400 1600 



m^rTT (MeV) 



bJO 




-400 



800 1000 1200 1400 1600 



m^rTT (MeV) 



Fig. 6cd 




Fig. 6ef 




800 1000 1200 1400 1600 

mTTTT (MeV) 




800 1000 1200 1400 1600 

mTTTT (MeV) 



Fig. 7ab 



bJO 




800 1000 1200 1400 1600 



m^rTT (MeV) 



bJO 




800 1000 1200 1400 1600 



m^rTT (MeV) 



Fig. 7cd 



bJO 




800 



1000 1200 1400 1600 



m^rTT (MeV) 



bJO 




800 



1000 1200 1400 1600 



m^rTT (MeV) 



Fig. 7ef 




Fig. Sab 




1600 



mTTTT (MeV) 



Fig. 8c 




Fig. 9ab 




1600 



mTTTT (MeV) 



Fig. 9c 



